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site moose model to describe the low energy Lagrangian of all same spin partner 
(LHT or UED like) models and check the gauge invariance of the amplitude. For 
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trast to those in supersymmetric models, there are spin correlations which affect the 
polar and azimuthal angle distributions of the quarks from the heavy partner U^^"^ 
decay. We show such effects would be visible in the i?Tmiss/^eff distribution and the 
reconstructed azimuthal angle correlation using MAOS reconstruction. 



Keywords: 



Collider Phenomenology, Spin Determination, Helicity Amplitude 



Contents 

|l|. Introduction |l| 

|2|. Quark partner production in a three-site moose model. |^ 

1^. Simplified Lagrangian and parton level distributions 

1^. Detecting spin correlations in jet level distributions [11 

|5|. Discussion and conclusion 16 



1. Introduction 

The existence of dark matter in our Universe is solid evidence of new physics beyond 
the standard model. Dark matter must consist of a new elementary particle, whose 
lifetime is much longer than the age of our Universe. One way to incorporate the 
dark matter candidate beyond the Standard Model (SM) is to introduce a parity 
structure. This parity structure can be incorporated in well-studied models such as 
supersymmetry (SUSY) (R-parity), little Higgs (LHT) (T-parity) JH, and universal 
extra dimensions (UED) (Kaluza Klein parity) p|, ^, ^. Under such a parity trans- 
formation, all the SM particles have even parity, their partners have odd parity, and 
the interaction vertices multiplicatively conserve the parity. The lightest parity odd 
particle becomes stable, and becomes invisible in collider searches. 

Observation of large missing transverse momentum at a collider would raise fur- 
ther questions about the properties of the underlying physics. The signatures of 
SUSY and same spin-partner models are quite similar 0, 0. In supersymmetric 
models, scalar quarks and scalar leptons have spin while in LHT and UED, quark 
and lepton partners have spin 1/2. The gauge boson partners are spin 1/2 in super- 
symmetric models, while they are heavy gauge bosons in LHT and UED. Therefore, 
it is critical to determine the spin of partners of the SM fields in order to distinguish 
different types of theories. To distinguish the spin structure of SM partners, we need 
to study a quantity that can be measured in the early stage of the Large Hadron 
Collider (LHC). At a hadron collider, the main production channels contain colored 
objects, whose main decay channels involve jets. The total cross section may serve as 
an initial hint to the spin of the odd partners ||^. However, at this stage of discovery. 
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small statistics make it impossible to reconstruct the mass, which the cross section 
depends on strongly. 

The spin-dependence of decay patterns of new particles is useful in determining 
the spin structure. The decay distributions depend on the spins of the initial, interme- 
diate and final state particles. For example, in the two leptons + jets + -Exmiss chan- 
nel at the LHC, the spin dependence appears in the jet invariant mass distributions. 
The simple mu distribution of dT/mu oc rriu arising from a g — )■ X2J ~^ XiJ^^ 
cascade decay in SUSY reflects the fact that the slepton is a scalar. In addition, the 
charge asymmetry in the rriji distribution arises from the polarization [|| of the 
spin 1/2 neutralinos. The polarization arises because the LHC is a proton-proton 
collider, so squarks are produced more often than anti squarks. 

The difference between the invariant mass distributions of SUSY signatures and 
those of same spin partner models have been studied previously [|10], |lT], [12], |13], |4 



The spin and interaction dependence is sizable, but such studies at the LHC may not 
be straightforward, because of the difficulty of selecting the correct jets and leptons. 
For a relatively clean jU signature, the decay branching ratio of g — )■ X2 is typically 
only around 30%, and X2 ^ ^ ~^ Xi is typically less than 10%. It is not the dominant 
decay channel and is therefore statistically limited. To distinguish the spin structure 
of quark partners in the relatively early running of LHC, we need to focus on the 
polarization dependence of the dominant decay channel. 

The other useful information lies in the spin correlation of particles from different 
decay chains. Here invariant mass distributions are not the useful quantity, rather 
the distributions of the polar angle and the azimuthal angle correlations |T^, |16|, ITT 



of decays. Due to the uncertainty of initial parton energy, it is not easy to see 
the angular correlation exactly at the LHC, but the effect should exist in various 
distributions in cases with strong angular correlations. 

In this paper, we study the spin correlations of the qq — )■ QQ process where 
Q is a spin 1/2 quark partner. We use a three-site deconstructed model which 
universally parametrizes low-energy effects in popular same spin partner (LHT or 
UED like) models with gauge boson and fermion partners. With such a simplified 
gauge symmetry breaking sector, it is straightforward to check the gauge invariance 
of the amplitude which guarantees good high energy behavior. In our model, the pair 
production and the decay of the spin 1/2 quark partner can have a spin correlation 
which is in contrast to the isotropic decays of scalar quarks. Among the possible 
signatures, the two high px jet + missing Et +X channel is promising, where X 
is an arbitrary number of jets and leptons. The quark partners can directly decay 
into a hght stable gauge boson partner and a quark, and the two high px jets come 
from their direct decays. It is essential that the jets coming from the two body decay 
tend to have a large pt because it makes their selection much easier. We stress that 
the channel we are proposing, which is one of the main discovery channels at the 
LHC, could be the most promising one to see spin correlations. Once the integrated 
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luminosity becomes closer to 1 fb^^, the LHC will be able to access partners with a 
production cross section above 0(1) pb~^. Therefore we propose simple distributions 
that are sensitive to the spin correlations, so that we would be able to determine the 
spin structure as early as possible. 

The paper is organized as follows. In Section 0, we introduce a three site model 
which gives us the low energy spectra of the massive gauge boson and fermion part- 
ners without model parameter constraints. We double check the gauge invariance of 
the amplitudes which guarantees their good high energy behavior. In Section we 
use two sample spectra in parton level simulations and present the polar angle and 
azimuthal angle dependences. We also calculate the amplitudes in different helicity 
channels. In Section |, we demonstrate how the jet level distributions are affected 
by the spin correlations. Section ^ contains our conclusions. 



2. Quark partner production in a three-site moose model. 

There are various models which contain SM "partners" with the same spin and an 
odd Z2 parity. Examples such as the universal extra dimension model ( "bosonic 
supersymmetry" ) [|], ^ ^, its variants with split mass spectra [|18|, |19|, little Higgs 



models with T-parities [|^ or even the warped models with KK parities belong to 
this group |2y]. In the process — )■ U\U\ followed by u\ decay into a quark and 
light stable particle (+ softer jets or leptons), one can safely neglect the effects 
of electroweak symmetry breaking, and all the massive particles gain their masses 
through the extra gauge symmetry breaking or the compactification with bulk gauge 
symmetry. From the dimensional deconstruction [^, ^ point of view, all of these 
models can be described by the three site moose diagram. Another reason we choose 
the three site moose digram here is because we need to describe a Goldstone boson 
with relatively simple mass mixings, which are needed in the i?^ gauge to check the 
gauge invariance of the amplitude and the high energy behavior. On the other hand, 
it would be very cumbersome if we were to use the the fifth component of the gauge 
boson in the language of extra dimensions. 

We consider a three site moose model which mimics the low energy theory 
for the scattering uu — )■ uxU\ in any kind of same spin partner (LHT or UED like) 
model. For the non-abelian case, our result would only differ by a color factor. The 
moose diagram is presented in Fig |I]. The gauge couplings in each moose are gB'-, 9a 
and gs-, and we set = 9b' so the model has a Z2 parity. The fermions charged under 
the gauge groups U{X)b and ^7(1)^' have equal Dirac mass terms —Mip^ip^ + h.c. 
and —Mijj^'ip^' + h.c. The two bifundamental scalar fields S and S', which we 
will call the "link" fields, are charged under gauge groups U{1)b', U{1)a and U{1)a, 
U{1)b with charge (1,-1) and (—1,1) respectively. They couple to the fermions 
through the Yukawa couplings with the same coupling strength yijj^'Eijj^ + h.c. and 
l/^/^fEV'l' +h.c.. 
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V'l' ^1 

Figure 1: The moose diagrams for the three site spht-UED model with a single 5D Dirac 
Fermion. An arrow into a site means that the particle transforms under the fundamental 
representation of the relevant site and an arrow out of a site means that the particle 
transforms under the anti-fundamental representation. The solid lines stand for Weyl 
fermions and the dashed lines represent scalars. The left panel shows how to get a left- 
handed chiral zero mode and the lower right shows how to get a right-handed chiral zero 
mode. 



The link fields get their vacuum expectation values (Sj^;) = (S^^) = u5jj. and 
spontaneously break U{1)b'i U{1)a and U{1)b into the diagonal group ^7(l)o- The 
kinetic term for the link fields Tr[(D^E)t(D/,E)] + Tr[(D^S')nD^S')] generates the 
mass terms for the massive gauge bosons. The mass matrix of the gauge bosons is 

{B^ B^) -gA9BU^ 2^>2 _^^^^„2 U J . (2.1) 
V -gAgBu' J \B'J 

We can find the mass eigenstates of the gauge bosons in terms of the gauge 
eigenstates 

Gl\ fs,/V2 c, s,/V2\ /BA 

^l] = \ -Vv^ 1/V2 . (2.2) 

GlJ \cJy/2 -s,cjV2j \B'J 

The inverse transformation between the gauge eigenstates and mass eigenstates is 

B,\ /sjV2-l/V2cjV2\ (Gl\ 

^ J = -s, \\gI \ , (2.3) 

B'J WV2 1/V2 cjV2j \GlJ 

where Sg = \^gA/g, Cg = gs/g and g = •\/2gi^ -|- g^. The corresponding gauge boson 
masses for G^, G^, G'i are 0, gsu and gu. 
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Similarly, the gauge invariant Dirac mass term and the Yukawa interactions give 
the fermion masses. The fermion mass matrix is 

/MlVlu \ /^|\ 
(^f ) Uif (2.4) 

V VLU mJ V^h 7 
We can also find the mass eigenstates of the fermions in terms of the gauge eigenstates 

= -1/^2 1/v^ , (2.5) 

\cf/V2 Sf CJ/V2J 

1 

-l/v^Ol/v/2 I Uil I . (2.6) 




I/V2 1/V2j 

The inverse transformation between the gauge eigenstates and mass eigenstates is 
>f\ fsf/V2-l/V2cf/V2^ 

^t] = [ -Cf 5/ I I I , (2.7) 









H 
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'0 -l/v/2 1/V2\ 

1 I I I . (2.^ 
1/^2 I/V2, 



Here, the fermion mixing angle is defined as Sf = V2yu/y/M^ + 2{yuy and c/ =M/ 
-\/iW^~+~2(yup, and the corresponding fermion masses for ip^, ip-^, •ip'^ are 0, M, 
+ 2{yuY. 

In this case, the model contains one gauge boson sector and one chiral fermion 
sector. The gauge boson sector contains one parity even massless state (0 mode), 
one parity odd massive state (1st mode) and one parity even massive state (2nd 
mode); The chiral fermion sector contains one parity even left-handed chiral fermion 
(0 mode), one lighter parity odd massive state (1st mode) and one lighter parity 
odd massive state (2nd mode). In order to mimic the SM and cancel the anomaly, 
one can also introduce another fermion sector which has a parity even right-handed 
chiral fermion as in the right panel of Fig. |l]. 

For the Goldstones, the interactions with fermions are 

D -iCfyijlip^-Ri + iCfyijj^iplTTi , (2.9) 
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where the symmetry breaking of a and a' is hnear 



S = (u + O") + ZTT . 

S' = (n + a') + in' (2.10) 

and eats the hnear combination tti = (tt' — n)/ 

Now we can derive the Feyman rules for the model and calculate the full ampli- 
tude. The relevant Feynman rules that we use in the calculation are presented in Fig 
||. The — — coupling is Qil = gCgSf/y/2, and the Goldstone tpR~ ~ 
coupling is g[]^ = —Cfy. For the scattering process uu — )■ uiUi, one must make sure 
that the calculated amplitude does not have bad high energy behavior. Indeed, as 
t'Hooft and Veltman in their proof of the renormalizability of the electroweak theory 
2^, it is the gauge invariance of the theory which guarantees the cancellation of bad 



high energy behavior of the amplitude ||2^. In order to check the gauge invariance of 
the amphtude, we adopt the general gauge and add the full amphtude including 
both massive gauge boson and Goldstone boson exchange. 

The amplitude coming from the t-channel exchange of the massive gauge boson 
in the gauge is 

iM = [-ig^] u{p')-i^PLu{p)— — ^— X (-) 

9'^ - - O) uik')^uPMk) . (2.11) 

We can extract the gauge independent piece from the gauge boson propagator as 

^fr-m.^(m. (2.12) 

So the additional gauge dependent piece can be written as 

iM^ = ]-{gBSffu{p')PLu{p) , ~\~'^} u{k')PLu{k) 
1 — ttiq^ rrtQ 

= {ycjfu{p')PLu{p) ^ ' u{k')PLu{k) . (2.13) 

Notice that mi = M = \/2yu/tf and md = gBU. 
The amplitude for the Goldstone interaction is 

iM^ = {-Cfyfu{p')PLu{p) 2 X (-) 

q TTT'G^ 

u{k')PLu{k) , (2.14) 

which cancels the contribution from the gauge dependent piece from the t-channel 
gauge boson exchange. 



- 6 - 




Figure 2: Feynman rules for the interactions used in calculating uu uiui. The mixed 
coupling strength for the gauge boson and fermions is giL = gCgSg/\/2 and the mixed 
coupling strength for the Goldstone and fermions is = —ycf. 



3. Simplified Lagrangian and part on level distributions 

With the above three site model, we can adjust the input parameters Ml, qa-, Qb, Vlu, 
Uru, Mr to get the free mass spectra for the massive gauge bosons and left-handed 
or right-handed chiral fermions and their massive vector partners. The relevant 
particles in same spin partner (LHT or UED like) models with the two jets plus 
missing energy signature are 

(1) a massive gluon partner Gh or electroweak boson partner Wh- 

(2) a light stable gauge boson Bh. 

(3) hght fermion partners Q^^\ Q^^\ L^^\ L^^\ (Here the superscript indicates the 

whether the fermion partner is partner of left-handed or the right-handed SM 
fermions, not the chirality of the partner.). 

All the fermion partner masses, gauge boson partner masses, and their interactions 
are free, but we assume the hghtest partner is Bh, the partner of the SM U{1)y 
gauge boson, and all partners have parity —1, so that Bh is stable. The other 
massive gauge boson partners, on the other hand, could be heavier or lighter than 
the massive fermion partners, depending on which model parameters we choose. The 
interactions involving gauge boson and fermion partners are 



-'int 



+giLB>tiYLQf^^^qL + giRB^nynQfim + (Lepton part) + /i.c] (3.1) 

Here gsi/B- 5'2Li and giL/R stand for the couplings between the SU{3)c, SU{2)l, 
U{1)y gauge boson partners and the SM left-handed or right-handed fermions and 
their partners. For simplicity, we take them to be the same as those of the SM gauge 
bosons in the following. We do not write the Goldstone interactions which ensure 
the gauge invariance of the amplitude. 
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Since the parton distribution functions of the first generation quarks are much 
harder than those of anti-quarks or gluons for protons, the production cross sec- 
tion of quark partners is much larger than that of antiquark partners. Gluon part- 
ners can decay into quark partners or anti-quark partners, however their produc- 
tion cross section is still subdominant. Therefore, light quark partner production 
pp Q^^'iQ^^yQ'^^^Q'^^^ is the dominant production channel at the LHC. Notice 
that the subdominant production of GhGh, GhQ^^^ and QQ may involve second 
level massive gauge boson exchanges in the s or t-channel unlike QQ production. 
The cross section for those subdominant process may be significantly different from 
that of our simplified model. 

The production and decay distributions of fermionic partners arc quite different 
from those of scalar quarks in supersymmetric models due to spin correlations. Here 
we demonstrate how those spin correlations appear in different helicity amplitudes. 
The amplitude for /i-) U^^^^{pf, hf)U^^'>''' {p'^, h'f) would be written 

as follows 



-9'' + ^ 

■X 



g2 - 



-|-gluon partner exchange contributions 

-|-cross diagram (3-2) 

where p's are momenta, /?.'s arc hclicitics of the particles, and a, h arc color indices. 
For the massive gauge boson propagator, we only write down the gauge independent 
piece since the rest of the pieces will be cancelled by the Goldstonc boson exchange. 
We only show the t-channel Bh exchange contribution and the g^^, and q^q" terms 
in the Bh propagator. The extension to the full amplitude is straightforward and 
for the numerical calculation we use the full amplitude. The amplitude arising from 
the Qfji^v term in the center of mass (CM) frame is 

^M, = ^9"YXMH^,^Jx„o J'^^, [25.„| /?/ - - Z?/)] , (3.3) 

2 q — 171^^ 1-12 J 

which is non-zero ior hi — h\ — 1/2 and Xf = hf — h'j: — 0. On the other hand, the 
q^q'^/m^ term after using the equation of motion becomes 

2 

in 

iM2 = -ig''^Y^Sab Sa'b'-r^ 2 ^h,{pf)PR UhiiPi) Uh'{p'f)PR Uhdp'i) , (3.4) 

which is enhanced by a factor of iTi^^Ry/rri^^. Note that this coupling is Yukawa- 
like, and it therefore flips the helicity of the fermions in the relativistic limit. The 
amplitude in the CM frame is 

2 / r- \ 2 

• /2-(A2t- r r r r '"( iflj / ^('M 

lM2 = -ig ^«^a6^a'6'5/,,,i5A,,o5Ai,0-^7-5— -2-^ 

'"'Bh\^ '"'Bh) 
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X 



2V,,-i/3/ + (-)'^^^(l -/?/)) 



(3.5) 



In the relativistic limit, (if — )■ 1, the amphtude is dominated by (/i/, h'j) = (1/2, 1/2) 
and (—1/2,-1/2) contributions. When mQ{R)/mB ^ 1, the h = —1/2 amphtude is 
dominant due to the Yukawa nature of the couphng. 

The decay of polarized Q^^^ (Q^^^ — )■ qBn) is not isotropic. The amplitude 
is dominated by the h = component of if rng^R) ^ mBu, therefore the decay 
distribution is essentially that of a spin 1/2 particle decaying into a light spin 1/2 par- 
ticle and a scalar. Because of the chiral nature of the vertices such as BHUR^y^PRUR, 
the quark in the final state is right handed. For a quark partner with h = 1/2, the 
decay amplitude is given by 



MviU^^^ -> uBh) oc e*'^e''Tr 



If^PR PH^Pr n (A + mum) 



where /i is the polarization vector with = —1 and fc// = —n ■ ks and Eb is the 
energy of Bh in the Q'-^^'s rest frame. Quarks from Q^^'^ decays with hg = 1/2 
tend to go in the direction of Q^^^ helicity as a result of helicity conservation. The 
distribution is proportional to 1 + cosOgg for hg = 1/2 in the massless limit, where 
9gq is the polar angle of the u momentum in the U^^'^ rest frame, where the z axis 
coincides with the U^^^ momentum in the CM frame. The decay distribution of Q^^^ 
{Q^^^ — )■ qBh or qWh) with h = —1/2 is the same as that of Q^^\ dV / dcosOgq oc 
1 + cos 9gq in the massless limit. On the other hand, it is opposite for an antiparticle 
decay. Because the number of produced quark partners is greater than the number 
of anti-quark partners at the LHC, the net polarization effect still remains in the 
signature. 

We calculate the amplitude using Madgraph/MadEvent assuming that the 
collision energy for the LHC is 7 TeV. The MAOS momenta parameter cards of 
the model are obtained by BRIDGE Fully spin correlated event distributions 
are generated by calculating the amplitude at the level of the final state partons 
qq — )■ QQ — )■ qqBnBn- Event distributions without spin correlations are generated 
for comparison by interfacing the qq — )■ QQ amplitude of Madgraph/Madevent to 
PYTHIA . We interface the events into AcerDet for detector simulation, but 
for jet reconstruction we adopt the Cambridge/ Aahen algorithm using FAST JET 
The number of events is O(IO^) before the cuts. For a typical production rate 



~ pb at 1 fb~^ integrated luminosity, the real number of events in the data sample 
would be about ten times smaller. Nevertheless, the qualitative features that we 
show in the plots would be the same. 
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In our simulation, we consider two sample points in the parameter space. We 
first choose the following spectra: The heavy gauge boson partner Gh and fermion 
partner Q^^^ have quasi-degenerate masses (we choose Q^^^ to be relatively lighter) 
and are much heavier than the lightest gauge boson partner Bh- In this case, one can 
easily identify the hard jets from a heavy Q^^^ decaying into Bh- We consider two 
different values of the light B^ mass, which has a large impact on the polarization of 
jets, as explained in the introduction, and illustrate how those different polarizations 
give us various different distributions. 

(A) mij(R) = 600 GeV, rrac^ = 700 GeV, tjibjj = 100 GeV. In this case, the sum 
of the production cross section of pp — )■ 

5.8 pb. The dominant production channel Q^^^Q^^^ has a cross section of 5 pb. 

(B) m^(R) = 600 GeV, = 700 GeV, but the lightest gauge boson mass is 
heavier, = 200 GeV. In this case, the production cross section is 2.7 pb, 
which shows that the Bh exchange contribution is important at Point A. 

For the spectra we choose here, the contribution of Gh exchange is larger than that 
of Bh exchange at point B, while the Bh and Gh exchange contributions are of the 
same order at the point A. 

In Fig. H, we show the cos6'Qq distribution for pp — i- U^^^U^^^ — )■ uuBhBh- We 
generate 10^ events using Madgraph, so each histogram contains 2 x 10^ entries. At 
Point A, the mass of Bh in the t-channel exchange is small, therefore the Yukawa- 
like coupling coming from the q^q^ /Tn\^ term in the Bh propagator is enhanced, so 
there are equal order (1/2, 1/2) and (-1/2, -1/2) helicity states for U^^^U^^\ which 
tends to decrease the overall polar angle dependence. In the Point B with larger Bh 
mass, the Yukawa coupling decreases. Only (1/2, 1/2) helicity states for JJ^^^U^^^ 
dominate and the overall polar angle dependence is stronger. 

We now discuss the azimuthal angle correlation in the uu — f/(^)f/(^) — ). 



uuBhBh distribution [|15], |T6|. The amplitude is expressed as the product of the 



production and decay amplitudes, 

iM{uu ^ [/(«)^7(^) ^ uuBhBh) = Y. Ph,h',-hMEcm, e)Dy,j,{e\ (P*)DMe'*, r) , 

h,h',h,h' 

(3.7) 

where Ph^h'JiJi' = ThyT^j^, and T^y is the helicity amplitude of uu — > U^^\h)U^^\h'). 
The decay matrix Di^j^{e*,(j)*) for U^^^ uBh is a function of 6* and 0*. Here 9* 
is the polar angle of the momentum of u in the rest frame of U^^\ where the z axis 
coincides with the momentum direction of U^^^ in the CM frame (the production 
frame), and 0* is the azimuthal angle relative to the production plane defined by 
the u and U^^'^ momentum directions. The decay matrix has simple azimuthal-angle 
dependence, 

D,,(r,r) = D,,(r)e^^('^-^)<^*. (3.8) 
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Figure 3: The cos Ogq distributions at Point A with = 100 GeV (left) and at Point 
B with ruBH = 200 GeV (right), tugh = 700 GeV and mg = 600 GeV. Ah histograms 
correspond to 10'^ events. 

The factor of exp{^{h — h)(t)*) is 1 if the production amphtude is dominated by a 
single hehcity state h = 1/2 or h = —1/2, but if both h = ±1/2 amphtudes are 
sizable, the azimuthal angle correlation is important due to the large interference 
term. In that case the distribution has the form oc 1 + acos(0) where is the 
difference of the two azimuthal angles. In Fig. ^, we show the distributions of at 
Point A [mBg = 100 GeV) and at Point B (m^^ = 200 GeV). The azimuthal angle 
correlation is stronger at Point A because the amplitudes for {h, h') = (1/2, 1/2) and 
{h,h') = (—1/2,-1/2) are roughly equal at this point. The polar and azimuthal 
angle dependences are complementary to each other. 

4. Detecting spin correlations in jet level distributions 

In this section, we discuss the effects of spin correlations in jet level distributions. 
Unlike at e~^e~ colliders, at hadron colliders it is not straightforward to see the effects 
by reconstructing the events exactly, because the collision energy of each event cannot 
be measured, and the lightest gauge boson partners are missed by the detector . In 
this section, we point out a few distributions which are sensitive to the effects. 

In Fig. H, we show the distributions of -ETmiss/^efF for the pp — )■ f/(-^)t/(^) process. 
In the left figure we show the distribution at point A; here the U^^^ decay is isotropic 
{^Bh = 100 GeV), which mimics SUSY squark production and decay. Fig. ^ (middle) 
is the fully spin correlated distribution for pp — > uuBuBh at Point A, and Fig. ^ 
(right) is the distribution at Point B (m^^ = 200 GeV). Here M^s is defined as 




(4.1) 
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Figure 4: The distributions of azimuthal angle difference (p of two U^^^ decays in the 
events at Point A (left) and Point B (right). All histograms correspond to 10^ events. 
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Figure 5: Distributions of ETmiss/Mes for pp — )• U^^^U^^^ — )• ururBhBh- Left; isotropic 
decays using PYTHIA at Point A. Middle; spin correlated decays at Point A. Right; spin 
correlated decay at Point B. All histograms correspond to 3 x 10^ events before the cuts. 

where the sum is over jets with pt > 50 GeV. Meg is correlated with the typical CM 
energy of the collisions. The missing transverse energy ETmiss is defined as 

^Tmiss = ^J {Pr'Y + {Pr')\ (4.2) 

where —Px^^) is the sum of the transverse momenta of the reconstructed objects. Up 
to the detector smearing and acceptance, it corresponds to the sum of the momenta 
of the invisible particles. 

The distribution of the events without spin correlation, (which is obtained by 
interfacing the pp — ?■ UU events to PYTHIA for isotropic decays) in Fig. |^ (left) is 
clearly different from the spin correlated decay. The spin correlated events tend to 
have smaller missing energy for fixed Mgg. The typical decay topologies at Points 
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Figure 6: Schematical description of event topologies of the process U^^^U^^^ — >• 
ururBhBh- Left: A case of polarized quark partners (Point B). The u quarks tend 
go in the same directions as their parent U^^^^s go. Right: the case with strong azimuthal 
angle correlation(Point A). Two u quarks tend to go in back to back directions in the 
production frame. 



A and B are shown in Fig. ^ (right) and Fig. ^ (left), respectively. The left figure 
shows a decay topology at Point B where the average polarization of the t/^^-* is 
right handed. The quarks tend to go in the directions of the parent t/^^^ momenta, 
therefore the energy of u (Bh) from U^^^ decay is boosted up (down) in average. 
The average i?Tmiss/^efr, which is correlated with PB^/mij(R), is smaller compared 
with that of isotropic decay. 

In Fig. 1^ (right), we show the other decay pattern reducing the average E'Tmiss/Miff, 
which is represented by the distribution at Point A. At this point, the average polar- 
ization of f/*^^-* is small, but the azimuthal angle correlation is stronger. In addition, 
the direction of U^^^ production is peaked forwards because mB„ /mjj(R) is small. 
Therefore the azimuthal angle distribution of u in the lab frame and that in the 
production frame are similar. The transverse momenta of the u quarks tend to be 
back to back in the production frame, therefore the two -Bj/'s tend to be back to back 
both in the production plane and in the lab transverse plane. The Bh transverse 
momenta cancel on average, so the total missing pj- becomes smaller. 

Next we discuss reconstruction of the azimuthal angle correlations. The simplest 
way to see this for our process is to look into the distribution of the leading two jets 
in the Lab frame. However, in this paper we use MAOS momentum reconstruction 

|31| to reconstruct the distribution. This is a method for estimating the direc- 



tions of missing particle momenta using the Mt2 parameter [^, which we explain 
below. 



The definition of Mt2 is expressed as follows p2 
Mt2 = iilin 



Pf+Pl" =PTmiss 



max ( mT(p2s,pf , "^test), mT(pS,p^, mtest) ) • (4.3) 



Here p^^ are the momenta of the two visible objects in the events. We adopt the 
two highest pt jets in each event, and pi^2 are those of the two invisible objects as 
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p^jg, and denote the calculated Mt2 by mT2 or mT2(2 jet). In this paper, we adopt 
minimal cuts to the jet rapidity rj and transverse momentum pt so that |?7j| < 2.5, 
Pti > lOOGeV, and pt2 > 50GeV. 

The minimization is taken over all pi and p2 that satisfy P1+P2 = Pmiss ^^nd 
{Pi? = {P2y = ml^,. 

We can define pill and pill from the (min) and p|'(min) which minimize the 
above equation. The conditions to solve piec are (In our case, Qi and Q2 are the 
intermediate massive heavy quark partners) 

p2 = („W I „(1)^2 _ 2 

(P^e^)' = ipielr = mLt, (4.4) 

for the transverse momentum 

PreF =Pr(min), p^eF ^^^(min). (4.5) 

For the events with mT2 ~ m™!^, pf (min) and p^(min) are very close to the true 
Bh transverse momenta. Therefore we define the reconstructed azimuthal angle 
difference 0(rec) based on Pq^ and Pq^, namely 0rec is the azimuthal angle difference 
in the CM frame of Pq^ +Pq2j where the z-axis is along the direction of the momenta 
Pq. There are two undetermined components in the prec and we have the two mass- 
shell equations for both sides of the decays, which give 2x2 = 4 solutions for 
(precjPrec). Two of them givc trivial solutions with 0(rec) = , and we study the 
other two non-trivial solutions, which lead to two azimuthal angle solutions, 0(rec, 1) 
and 0(rec, 2). 

The merit of using MAOS reconstruction is limited. In Fig. we plot the differ- 
ence between the reconstructed 0(rec, i) of the parton level event and the true parton 
level azimuthal angle difference (/)(true) for the events near the end point with 771^2 > 
500 GeV. The plot is based on 20000 generated events where ~ 2500 passed through 
the cut. Here the solid histogram is A0(best) = minj=i^2 |(0(rec,z) — (/)(true))|, the 
dark dotted histogram contains both of A0(i) = |0(rec, i) — 0(true) | (scaled by a fac- 
tor of 0.5), and the thin dotted histogram is A0(Lab) = ||0i — 02| — 0(true)|, where 0i 
is the azimuthal angle in the Lab frame. The matching between 0(true) and 0(best) 
is good. However, once we take into account both of the solutions, 0(rec, 1) and 
0(rec,2), the improvement from the simple Lab frame quantity A0(Lab) is moder- 
ate. The mean value of A0(best) is 0.28, while it is 0.41 for A</)(rec, i) and 0.50 for 
A0(Lab) respectively. 

At jet level, we have to consider the possibility that the selected momenta do not 
match with parton level momenta. In Fig. |](left) and Fig. |]( middle), we compare 
parton level and jet level distributions in 0(rec)-0(true) plain. In the left plot, we 
plot 0(rec) with 0(true) at parton level, and in the middle plot (/)(rec) is calculated for 
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phi(best)— phi(true) 



Figure 7: The difference between reconstructed azimutlial angle difference and (/>(true) at 
parton level. The solid line shows the difference from the best solution, the thick dashed 
histogram plots both of the reconstructed solutions, and the thin dotted line is the difference 
to the lab frame azimuthal angle. All histograms corresponds to 20000 events before the 
cuts. 



the two highest pt jets in the events after the detector smearing. For 0(true) ~ the 
probability to mis-reconstruct the azimuthal angle is high. This is because when two 
partons go collinear, the correlation between parton momenta and jet momenta is 
not good. In Fig ^, we compare the distribution of 0(rec) at Point A (solid line) with 
the distribution of the events without spin correlation (a thick dashed histogram) 
and the parton level distribution without spin correlation (a thin dashed histogram). 
All distributions have a peak at = tt, but the distribution with spin correlation is 
most strongly peaked at </)(rec) ~ vr as expected. 

The distribution of 0(rec) is significantly altered by the kinematical cuts to 
reduce the background. To explain the effect of the cuts, we show the correlation 
between 0(true) and mT2 in Fig. ^j. In Fig. |^ left (middle) we show the distribution 
with (without) spin correlation. Here, mT2 is calculated for the two highest pt jets as 
the visible objects of the events, and mtest = f^BH- We plot both of the reconstructed 
azimuthal angle solutions in the plot. The expected endpoint of the mT2 distribution 
is 600 GeV. For the two jets + missing signature, it is known that requiring a lower 
limit on the value is an efficient cut to remove the standard model background 
Q. In the region ~ 0(7r), m,T2 tends to be large (small). At Point A the number of 



events at (prec ~ is small compared with the distribution without spin correlation. 

In Fig. I (right), we show the 0(rec) distribution with various mT2 cuts. Solid 
and dotted lines with a sharp cutoff around 0(rec) < 2 correspond to mT2(2 jet) > 
400 GeV, where mT2(2 jet) is the 771^2 calculated with the two highest pt jets as the 
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phi(rec) phi phi(rec) 



Figure 8: Left: parton level (p{iec) vs. 0(true) distribution at Point A. Middle; jet 
level (p{iec) vs. (p{true) distribution at Point A. Right (j){rec) distributions at jet level 
(solid histogram), jet level without spin correlation (thick dashed histogram), parton level 
without spin correlation (thin dashed histogram). See the explanation in the text for 
details. Each plot corresponds to 2 x 10^ events before the cuts. 



visible objects. About 3600 events remain after the cut for spin correlated events 
and 5000 events for spin non-correlated events. Significant suppression is observed 
for (/)(rec) < 2 for the events with spin correlation. The sharp edge at ~ 2 appears 
because the events above 0(true) > 2 for 771x2 > 400 GeV are kinematically forbidden. 

The other two histograms with significant tails beyond 0(rec) > 2 are for the 
events with mT2 >400 GeV, but the mT2 is calculated for all jets in the events using 
inclusive mT2 defined in [^, 35 1. Here the jets and leptons in an event are clustered 



into two visible objects to obtain mT2- The number of events remaining after the cut 
is ~ 5300 for spin correlated events and ~ 6500 for spin non-correlated events. The 
tail events are the QQ events with significant initial state radiation. The events with 
771^2(2 jet) < 400 GeV contaminate the 0(rec) > 2 region due to the initial state 
radiation. 

The number of events with </)(rec) > 2 is roughly the same between spin corre- 
lated and non-correlated samples, because ISR and the leading jets from Q decay 
are not spin correlated. This is encouraging because the events in this region may 
be regarded as control samples to determine the total number of produced events. 
Together with the number of events with 0(rec) < 2, we may be able to determine 
how much azimuthal angle correlation exists in the events. To judge the usefulness 
of the method, it is necessary to conduct a more careful study of the dependence on 
more general production and decay patterns, and the effects of other cuts such as 
missing Et- 

5. Discussion and conclusion 

In the coming years, the LHC will accumulate more events, and the signature of 
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Figure 9: The (/)(rec) vs mT2 distributions with and without spin correlation (left and 
middle plots, respectively) Right: The 0(rec) distribution with spin correlation (solid his- 
togram) and without spin correlation (dashed histogram). Each plots generates 2 x 10^ 
events before the cuts. 

supersymmetry, namely an excess of events with may be found. However, there 
are other models with similar signatures but different spin structures. In this paper 
we discuss a way to measure the spin dependence of SUSY-like signatures, especially 
those that might be observed in the early stages of the ATLAS and CMS experiments. 

We focus on the spin correlation in the 2 high px jets + missing Et +X signature 
at the LHC. The signature would arise from a pp QQ process when 1) the model 
has a colored quark partner Q that directly couples to quarks. 2) The colored 
particle dominantly decays into a lighter parity odd particle x ^-nd a quark q. The 
X are gauginos or Higgsinos in SUSY or heavy gauge partners in same spin partner 
(LHT or UED like) models. 3) The mass satisfies tuq 3> m~^. The x ^lay further 
decay into other particles involving the lightest "parity odd" particles. In that case, 
the jets from Q decays are prominent among the jets and initial state radiation 
Depending on the spin of Q and its polarization, the two jets show non-trivial 
spin correlations, both in the polar angle dependence of Q qx decay and in the 
azimuthal angle correlation of the two quarks. The effect can be seen in -Exmiss vs 
Mgfj distributions, and the reconstructed azimuthal angle correlation using MAOS 
momentum. 

The expected spin correlation depends strongly on the production process. In 
this paper we have studied the production process based on the t channel exchange 
of the heavy gauge bosons expected in same spin partner (LHT or UED like) mod- 
els. In those models a parity odd quark Q couples directly to a quark q, and the 
interaction among Q, q and parity odd heavy gauge bosons Bh and Gh could be 
chiral. We have shown that the helicity structure of the amplitude is very sensitive 
to the mass and the interaction of the heavy gauge bosons. The enhanced helicity 
h = component plays an especially important role. As an example, we show the 
distributions of uu — ?■ U^^^U^^\ The polar angle dependence comes from overall 
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polarization, while azimuthal angle correlation comes from interference between spin 
flipping Bh exchange and spin conserving heavier Gh exchange. The azimuthal an- 
gle correlation is visible because the amplitude is forwardly peaked due to the Bu 
exchange contribution. 

Due to the azimuthal angle correlation of the events, the number of events may 
be suppressed near the (2 jet) end points. The high region is the signal 
region because the SM background is smaller, therefore spin correlations affect the 
estimation of the total production cross section. We note that there are phase space 
boundaries of the events of pp — )■ QQ —> uu + -ETmiss- Some events can lie outside 
the phase space due to ISR, and they are less affected by the spin correlations. We 
note that such events may be actively used to estimate the total cross section of the 
events. 

In this paper, we do not study the corresponding distribution in SUSY models 
in detail. In supersymmetric models we do not expect spin correlations because 
the partner of the quark has spin 0. This leads to a high i?Tmiss/^efr, a hard 
end point, and a flatter reconstructed azimuthal angle correlation for two jets. Of 
course there are other important differences, the squark squark pair production is 
generally small in SUSY-like models, and gluino squark co-production is dominant. 
The gluinos may decay into squark and quark, and the ISR of the process involving 
the gluino is larger, so selecting the correct jets coming from squark decay may not 
be straightforward. However, it is shown in Ref. [3C] that the two highest pt jets in 
the events are likely from q XQ decay. The inclusive distributions with high 
Pt jets show a sharp end point for the model parameter rriq < rrig. This means the 
two jets from squark decays are the dominant part of the event activity. Therefore 
we believe the feature discussed in this paper would be useful for SUSY studies as 
well. 

For extracting the physics behind the signatures at the LHC, using fully spin 
correlated amplitudes for production and decay is essential. To reproduce the struc- 
ture of the amplitude, the amplitude of qq — )■ QQ — >■ qqBuBu + X should be fully 
calculated. In this paper we calculated the full amplitude using Madgraph p5[, but 
other amplitude calculators or generators such as Comphep/Calchep [0 and Her- 
wig++ |]SS| should be able to reproduce the effect. On the other hand, when we 
interface the two to two process to PYTHIA the amplitude does not have the 
proper spin correlation, and therefore it is not appropriate for this study. 
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